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Abstract 

In Einstein theory of gravity the initial configuration of metric field and its time 
Q\ ' derivative are related to matter configuration by four equations called constraints. 

We use the method of conformal metrics (York Method) to solve constraints and 
O , find an analytic set of consistent initial data for linearized Einstein field equations 

in a perturbed constant curvature space-time. They are explicitly covariant and 
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more compact than decomposition of quantities to scalar, vector and tensor. This 
method is independent of type and physics of matter fields and is extendable to 
higher-order perturbative calculations. As an application example, we apply this 
method to two commonly used matching conditions during a phase transition and 
compare and interpret the results. 
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By a set of initial conditions, we mean a configuration for matter, radiation, and metric 
fields, and their time derivatives on a given initial 3-space. This information is neces- 
sary for solving Einstein equations. The task of defining this configuration is not trivial 
[Lichnerowicz 1944]] , ||Arnowitt et al. 1962|| . In contrast to Newtonian theory, due to dif- 
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feomorphism gauge symmetry, the initial value of the metric and its time derivative can 
not be specified in an arbitrary way. The configuration of matter fields partially deter- 
mines them. 

To define initial conditions in general relativity properly, one must define a time slic- 
ing of the space-time i.e. a diffeomorphism between the space-time and a 3+1 manifold 
York 197211 , ||0 Murchadha fc Yorkl97l , flChoquet-Bruhat fc York"l98C|l , jDurrer 1994 



From this operation, in addition to evolution equations for geometrical quantities, one 
obtains the restriction of the Einstein equations on the space-like 3-space component of 
the 3+1 manifold. These equations don't evolve with time and for this reason they are 
called constraints [|Arnowitt et al. 1962j| . The initial condition for geometric and matter 



fields is defined as their configuration on a space-like 3-space. As constraints are not 
dynamical equations, it is necessary and sufficient that field configuration on the initial 
3-space satisfy them. 

In cosmology, specially for studying the evolution of small perturbations, it is customary 
to decompose Einstein equations, as well as energy-momentum tensor, to scalar, vector, 
and tensor components. Evolution equations and constraints for each type of fluctua- 
tions are solved separately. This procedure drops some of components from constraints 
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and makes them easier to solve, specially when only one of components, usually scalar 
component, is studied. Nevertheless, evolution equations and constraints remain coupled 
and one has to solve them together. Moreover, the resulting equations are not usually 
explicitly covariant (although it is possible to perform the decomposition in a covariant 
manner [Stewart 1990|| ). 

Having an exact solution of constraints is also important in numerical solution of evolu- 
tion equations. Usually, it is very difficult to keep the conservation of energy-momentum 
in a numerical calculation. The exact solution of constraints assures that the initial data 
satisfies the conservation laws. It also can be used at each step of calculation for check- 
ing/correction of conservation violation. 

The mathematical aspects of the initial conditions for the Einstein field equations have 
been studied and clarified in an outstanding work by J.W. York on the initial condition 
problem and its relation to conformal gravity ||York 1972| , ||Q Murchadha fc York 1 971 1. 
This formalism separates initially constrained components of metric and extrinsic curva- 
ture, and allows a more detailed insight to the physical nature of unconstrained compo- 
nents without knowing anything about the matter content of the theory. 
Here we apply this method to linearized Einstein equations, solve constraints analytically, 
and find a consistent set of initial conditions for perturbed constant curvature cosmogo- 
nies. The results are independent of details of the physical model under consideration. 
As an application example, we use the results of this method to discuss matching during 
a phase transition. 

We first briefly remind the mathematical formulation of the problem. Then, we explain 
analytical solution of constraints for flat and constant curvature perturbed space-times 
and finally we apply the solution of constraints to matching conditions on a phase tran- 
sition surface and we obtain the initial values allowed for unconstrained components. 



2 3+1 Gravity 

We assume that by a diffeomorphism transformation, the space-time is divisible to a 
space-like 3-manifold and a one-dimensional time like space (curve). The general form of 
the metric is: [] 

ds 2 = -a 2 drj 2 + g ij (dx i + p i dr ] )(dx j + (3 j dr]) dt 2 = a 2 drf (1) 

7] is the conformal time. a(j],x) and f3(r],x) are called lapse function and shift vector (in 
3-space). The 3-tensor is the induced metric of the 3-space. 

Einstein field equations can be expressed by a set of first-order differential equations 
which depend on the field set {a, [3\ gij, Kij}. The extrinsic curvature is defined 
as the covariant derivative of the normal vector to the 3-space = —n^ = — a^T°- 
Misner et al. 1973"|. The vector n M is the normal unit vector of the 3-space. With respect 



1 Greek indexes correspond to space-time and Latin indexes to 3-space. We use units where c = 1 and 
% = 1. The notations "|", and ";" are used respectively for partial derivative, covariant derivative in 
3-space and covariant derivative in space-time. 
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to these fields, Einstein field equations take the following form |[Piran 198C|| : 



dr] 
dr] 



-2ag ik K k + (3 k g ijM + g ik (3 k + g jk (3% (2) 

p k Kl k + K k ^ k ^ - K k (3\ k - oP\i + 
a{KKi + ®bJ + 8nGgi k T lk ). (3) 

{3) R- Ki j K/ + K 2 = 16ttGT*. (4) 
K^-K\i = SirGT*. (5) 

The quantities T* = T^n^B" and T* = T^ u n^n v can be interpreted respectively as 
projection of energy-momentum flux on the 3-space and on its normal. B l is the projection 
operator on the 3-space, B^n^ 1 = 0. For metric ([]]): 

n» = oT x {\,F). (6) 
T* = a-\T 0l -T lk [3 k ). (7) 
T; = a- 2 (T 00 -2T ^ + r^W). (8) 

Equations @) and @ are dynamical equations for the evolution of the 3-space. Equa- 
tions (Q) and (^) don't have explicit time dependence. They are constraints. If they are 
satisfied by the metric and extrinsic curvature at the initial time to on the 3-space, they 
stay valid for ever. 

Gauge symmetry allows to fix arbitrarily the value of at(t,x) and (3(t,x). In synchronous 
gauge that we will use in this letter, a depends only on t and (3 l (t,x) = 0. The gauge 
symmetry assures that these relations remain valid after the evolution of the dynamical 
system. In fact, 4 of 16 equations (0) - fl5|) are automatically satisfied due to Biancci 
identities. This results to an equal number of equations and fields. The physical reason 
behind existence of constraint equations is the gauge symmetry. A system with gravita- 
tional interaction is a constrained dynamical system and the initial value of g^ and 
can not be arbitrarily chosen. The initial g^ and fields configuration on the initial 
3-space depends on the matter configuration and must satisfy constraints @ and (||). 



3 Conformal Metric Method 



In || York 197T| and [ p Murchadha fc Yorkl974]j it has been proved that unconstrained de- 
grees of freedom of a gravitational system on a space-like 3-space are conformally equiv- 
alent. This allows to separate constrained and unconstrained components of gij and Kij. 
The extrinsic curvature tensor of the 3-space can be decomposed to traceless-transverse, 
longitudinal (traceless), and trace components: 



K, 



ij 



S {J + (LW)ij + -Kgij. 



(LW) 



Kl 



w Aj + w. 



J I* 



(9) 
(10) 



The scalar K is the trace of and 5^- is its traceless- transverse component (Sfj = 0). 
W l is a 3- vector field that generates the longitudinal (traceless) part of Kij. 
The conformal transformation of a 3-metric gy is defined as gij = (j) A gij. We call gij 
the base metric and the physical metric which satisfies Einstein equations. Under this 
transformation, @ and (0) keep their form. The only difference between base or physical 
extrinsic curvature decomposition is the use of base or physical metric in the contraction 
of indexes. Therefore: 



K tJ = § ij + (LW) ij + -Kg ij . (11) 
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{LW) l3 = W i{j + - \giiW^. (12) 

t) k = r} fe + 20- 1 (^| fe + 4^--^). (is) 

The trace K and the vector potential W % are considered as invariant under this transfor- 
mation. The traceless-transverse component Sij is chosen to transform as: 

Sij = r 2 Sir (14) 

(When only one initial 3-space is considered, the conformal transformation of is ar- 
bitrary. However, it has been recently shown that by assuming two infinitesimally close 
3-spaces and a mapping between them, this transformation rule is imposed by Einstein 



equations [[York 1998fl ). Even though the choice for transformation of T* % and T* is not 



unique, following choices guarantee the existence of a non-spacelike energy flow from the 



3-space ||Q Murchadha fc Yorkl97l : 



T* = 0- lo T". (15) 

?; = r s Ti (i6) 

With these relations, the longitudinal component must transform as: 

{LWf = <f>- 4 {LW) ij . (17) 

cjij and Kij satisfy the constraint equations (|J) and (|5|) if W l and 4> satisfy following 
equations: 



(Lwy% 




(18) 


-8A0 


= -R<f> + Mtt^" 1 + 2M TL (j)~ 1 + 






(M L - ^K)^ + IGttGT;^ 3 . 


(19) 


M TT 


= 9ac9bdS ab S cd . 


(20) 


M TL 


= gacg bd S a \LW) cd . 


(21) 


M L 


= g ac g bd (LWy b (LW) cd . 


(22) 



All indexes and derivatives are contracted with the base metric g^. Equations (]T8[) and 



(19) are obtained from constraint equations (HI) and (EI). These equations show that 4> and 



W % are the real constrained degrees of freedom of the initial configuration of fields. The 
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unconstrained degrees of freedom are gij, S{j and K and they are all conformally related 
to their physical counterpart. Note also that this formulation of the initial value problem 
is gauge independent. 

The real difficulty in solving equations flTSP and ( |TP| ) is that they are highly nonlinear 
and in general completely coupled. This fact has encouraged the use of other methods, 
specially for nonlinear numerical calculations ||Piran 1980 . 



4 Initial Data for Linearized Einstein Equations in 
Flat Cosmogony 

We consider an initial equal-time 3-space S in a flat universe with small perturbations 
(it is always possible to redefine coordinates such that the 3-space become equal-time). 
A priori S can be any hypersurface, but interesting cases are those with constant K or 
K = Khomo + SKj such that Ku = (SK)\j and Khomo is regarded as the trace of for a 
3-space obtained from mapping S to the homogeneous background manifold. In fact, K 
is an appropriate quantity to time label a spacelike 3-space because it has a non-negative 
derivative in the direction of timelike orthogonal vector to the 3-space fYork 1972fl . Other 



physically interesting 3-spaces e.g. one with constant energy in a perturbative theory are 
close to a constant K 3-space. Here we only consider these types of 3-spaces for imposing 
initial conditions. 

We take the base metric to be the metric in a flat homogeneous Freedman-Lemaitre 
cosmology, i.e.: 

= a 2 (t)5ij. (23) 

This choice is physically motivated, because for a homogeneous universe, this metric 
satisfies Einstein equations. In addition, it simplifies all calculations. 
With this choice, in a flat universe with small perturbations, the conformal factor will be 
close to 1: 

0(x,t) = l + <ty(x,*) (24) 
For linearized Einstein equations in synchronous gauge, the physical metric is: 

g^ = a 2 {t)(5ij + hij) 

= a 2 (t)0 4 (x, t)8ij = a 2 {t)5 i3 {\ + 450(x, t)). (25) 

= 4Sij6<f>. (26) 

This special form of metric fluctuation is not a restriction of perturbations to scalars and 
is not equivalent to Newtonian gauge because its application is uniquely on the initial 
3-space. As we have seen in the previous section, the transverse-traceless component of 
the extrinsic curvature is an unconstrained quantity. If on the initial 3-space it is not 
zero, this pure tensorial curvature perturbation will contribute to g^ evolution (see (0)) 
and induces a purely tensorial metric perturbation. 

Another important point is that it is always possible to choose a coordinate system on 
the 3-space with a metric like (|25|). It is well known that synchronous gauge does not 



define the gauge completely and allows a redefinition of the coordinates which preserves 
the gauge. In a gauge preserving coordinate transformation like the following: 



t' = t and x* = x l + aV + £ ij Uj (27) 

the arbitrary fields ip (local translation) and u>i (local rotation) can be chosen such that 
§ij gets the form given in (p5|). After calculating the initial condition, one can return to 
the previous coordinates or to any other gauge. 

From definition of projected energy- momentum tensor Equation. (§), in synchronous gauge 
T" = Tq. In an approximately homogeneous and flat universe with small perturbations, 
Tq = 0(1) in the comoving frame. Therefore, T* % = + ST* 1 . Regarding (|T5|), at 
first-order of approximation, T* 1 = T*\ From @, one can conclude that at zero-order 
Kl- = 0. This means that at zero-order, the physical extrinsic curvature is transverse, 
and therefore, in spaces with small perturbations, the value of its longitudinal part must 
be small. Equation ( |TTD shows that the longitudinal part of the base extrinsic curvature 
also must be small. 

For small perturbations, equation flTE|) becomes: 



[(1 + 65(f)) {LW)% = -8nGT* 1 + -(1 + 68(f))K^. (28) 

3 

The term 5(j)(LW)^ is of second order and one expects that for smooth fluctuations, its 

derivative must be negligible with respect to (LW) 1 ^. As mentioned earlier, we assume 
that = 0(1). Therefore, if: 

mLwy% « (LW)\h ( 29 ) 

(ESI) reduces to: 



{LW)% = -8nGT* 1 + — ^ SnGT** + — . (30) 
IJ 3 3 

The solution of this equation in the case of a flat base metric is trivial. With a variable 

change x n = a{t)x\ (PD| ) becomes: 

d' j d'W l + Id^d'Wl = -87rGT*Xx', t) + (31) 
3 3 

d' means partial derivative with respect to x'. The solution in real space is: 

2nG 



(2n) 3 J k 4 
h l 5K 

[k 2 (4T»(k,t) - —) - k%T*i(k,t)\. (32) 

For solving equation (0), one needs also Sij the transverse part of the base extrinsic 
curvature. It is an unconstrained component and must be chosen according to the physics 
of the system. To understand better the physical role of SV, , note that its definition: 

Sg = (33) 



can have non-trivial solutions. It has been proven that this equation allows a solution 
with a singularity called monopole solution (a black hole, a point like object of finite mass, 
or a geometrical singularity at origin) Powen fc York 1980|| : 



-{e mn S m u n u> + e^ mn S m u n u l ). (34) 



The vector r is the radial vector through origin and u l is the radial unit vector, P % and 
S l are respectively the linear and angular 3-velocity of the singular point. This solution 
can be extended to a solution with iV sine ularity at fW, % = 1 • • • N ||Bowen E York 198(J|| 



Brandt fc Bruegmann 1997 1 . In fact, this is a solution for ([30]) in vacuum. However, in 
this case, (LW) y ' is transverse and it is possible to add ([34]) to any solution of (|33|) or in 
general to K^. The existence of this solution reflects the effect of a velocity field on the 
extrinsic curvature and on the evolution of the metric (equation @). Note also that T^ v 
does not fix the velocity field. This make a solution of type ([54]) independent of fl3"2|). In 
decomposition method, the linear and angular velocity fields appear in the decomposition 
of Tjjj, and their geometrical role is not as evident as here. 



The solution (|34"D can be generalized to a non-singular matter distribution: 

3 



J V— \x\ 



v-{x} v a 2\r — r'\ 2 
(pVK + p*(rV - (g ij - uV> c (r'H) + 

v-{x} \r — r | d 

(e im " Sm (r'K^' + e^ Sm (r')^). (35) 

The integration is performed in the horizon of point x except the point itself. Vector u % 
is in the direction of the source at r*. The vectors p l (r) and s*(r) are linear and angular 
velocity field densities with following definitions: 

e^% k = 0, 4 = 0. (36) 

In general, these vectors can be obtained from matter distribution. Boltzmann equation 
relates them to non-gravitational interactions [ JEhlers 1971| ], 



Finally, S ZJ can have a purely transverse-traceless component independent of matter dis- 
tribution and related to relic gravitational waves. Therefore S 1 ^ = S l J el + S l r 3 elic . 
Having S l i for a given distribution of matter and relic perturbations, we can now use 
equations ( PD| ) to ( P^P to determine the coefficients of (|T9"D. Using approximation 



and the fact that for the chosen base metric = 0, equation ( |i~9"| ) changes to: 
A5cf) = -~[(M TT + 2M TL + M L - 

o 

2 -Kl omo + lQnGT: homo ) - (7M TT + 2M TL - 5M L + 
10 4 

—K 2 + AMGT: homo )]5<t> + IGttST: - -K homo SK. (37) 
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In a homogeneous universe: 



9ij = a 2 5 {j and K i:j = -^-<%- (38) 

From (|4]), one can see that: 

2 2 # 2 2 

\QnGT^ homo — -K homo = lQixGT^ homo — -K homo = 0. (39) 

To keep 5<p small according to our small perturbation assumption, the matter has to 
have a small average momentum. This makes all coefficients Mtt, Mtl and Ml of 
second-order and negligible. It is also another demonstration of power of this method. 
|Veeraraghavan fc Stebbins 1990|| also arrives at the same conclusion for the vorticity of 



a perfect fluid. Here this result is obtained without any assumption about type and 
state equation of matter. At first-order, the velocity perturbations contribute only to the 
evolution of the metric (equation (Q)) and not to constraints. 
With these perturbative simplifications, (|37|) reduces to: 



AS(f> = -2nG5T: + -K homo 5K = (-K 2 + 6ttGT; Wo )^ - 2nG5T: + -K homo 6K. (40) 

and can be solved as: 

+ , 1 r , 3 ,2nG5T:(x',t)-lK horno 5K(x>,t) 

1 f « ]r - ta( t)s,,k^ ^G5T:(k,t) - lK homo SK(k,t) 

l z7r J J K ^ 12 JX *homo^ un '^ 1 *horno 

(r and r 1 are as in (|35D ). From the form of the base metric, T* homo = T* homo . 

The first expression for 8cj) shows that physical metric gij is similar to the metric for scalar 

perturbations alone (see e.g. ||Bertschinger 199l| ). It is the effect of special coordinate 



choice on the 3-space. 

We now have all quantities necessary for definition of a set of initial data and solution 
of evolution equations (0) and ([J). The essential characteristic of linearized constraints 
is that equations for W l and (p are decoupled. In fact, for all order of perturbations, 
Sep decouples from (ESI) and therefore, this method is easily applicable to higher-order 



perturbative calculations. 



5 Non-Flat Cosmogonies 

For the general case of a space-time with constant curvature, we can use a flat base metric 
as before. The perturbative expansion of (p will take the following form: 

0(x,t) = o (l + ^(x,t)). (42) 

0o = i r- ( 43 ) 

(1 + fr^ 
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r 2 = 5 ij x i x j (44) 
(jij = « 2 (t)(7 ij + hij) 

= a 2 (t)0 4 (x, t)5 l3 = a\i) ll3 {l + 4<5</>(x, t)) (45) 

(1 + T r 2 ) 2 

K is the constant curvature of the homogeneous universe. After performing a variable 
change from x to x' as before, the equation for W l will become: 

d fi W j + d' j W i - -S ij d'W c d H K 

(l + fr'2)3 J 3 1 ; 

where A"' = To solve this equation, we assume, without loss of generality, that W 1 

can be decomposed to: 

W i = ipV^^U. (48) 
V = (1 + ^r' 2 ) 3 . (49) 



Putting ( [4lf) into (|47|) gives a system of equations for V % and [/: 

+ d' ] d' 3 v i - ^d'fty = -sttGt^x',*) + (50) 



3 v 7 y 3 

3' 4 [7 = (51) 



F ij = V s + V*#ty + h ij V%ijj. (52) 

o 

Equation fl5"0|) is exactly the same as equation ([H|) and therefore its solution is the same 
as fl32"|). In (0), the right hand side is known once fl50D is solved and the solution of (|5T|) 
is: 

"(M) = — / # fe ^ W «-3M|^W (53) 

Equation fl3*T|) changes to: 

A50 = (^K 2 homo + 6nGT: hom MH - 2nGfa 3 5T; - ^ 5 K homo 5K. (54) 

To solve this equation, one can expand 5(f) and ST* with respect to Spherical Harmonic 
functions: 

5(f) = ^a m i(r)Y m i(9,ip) 5T* = ^T m i(r)Y m i(6 , ip) 

ml ml 

5K = ^K ml (r)Y ml (e,<p) (55) 

ml 

This results to the following equation for a m i. 

r 2^L +2r ^i + W + 1) _ r 2 Air))ami = 



dr 2 dr 

— [r 

dr dr 



'' ( 'A + (l(I + l)-r 2 i(r))fl m! = r 2 (B(r)T ml + D(r)K ml ). (56) 



o 



A(r) = « 2 (t)(^L mo + 6^GT; /lomo )^ (57) 
B(r) = -2nGa 2 (t)(j)o 3 (58) 
D(r) = ~4K homo . (59) 

This equation has a polynomial homogeneous solution. For simplifying our notation, in 
the following we ignore m and I indexes. We assume: 

a(r) = Y / d n r n . (60) 



for the homogeneous solution. Replacing a m i in Q5"fip with this expansion, one obtains the 
following recurrent expression for the coefficients d n : 

8 n 

dn = Y, 7?[(n - 2% + 2)(n - 2* + 3) + 1(1 + l)]4-2i +2 . (61) 
i=o Gl 

Ci = « 2 (t)(^^ omo + 6vrGT; Wo ). (62) 

C 2j = - (63) 
i!(8-i)! v ; 

One can construct two independent solutions for two boundary conditions: 

1) \a\ < oo for r -> rf n = 0, n < 0. (64) 

2) |a| < oo for r — > oo =>- d n = 0, n > —16. (65) 



As the derivative terms of (p6|) are complete, a Green function can be found for this 
equation ( ||Zwillinger 1989(1 ) • We ca U the above homogeneous solutions a^(r) and a^(r). 



The Green function and the complete solution of a m \ are determined as: 
G[r-r') = 



° (1 W 2 V) <r<r', 



r' 2 w(r') 
a (l)( r ') (2)( r ) , 

U M = a (1) (r)— a {2 \r)——. (67) 

ar ar 

a m i{r) = J dr'r' 2 (B(r')T ml + D(r')K ml )G(r;r'). (68) 
This completes the solution of constraints for perturbed constant curvature spaces. 



6 Matching Condition 

The results obtained in previous sections can be applied to any perturbative cosmological 
context. An straightforward and interesting application is the determination of initial 
metric and extrinsic curvature perturbations in a universe filled with species which have 



1 n 



a distribution /(x, t, p) in their classical phase space. If they have mutual interactions, 
one usually has to solve the Einstein-Boltzmann equations numerically. Division of these 
equations to scalar, vector and tensor components increases significantly the amount of 
equations to be solved and consequently the calculation time. This issue and the detail 
of calculation will be discussed elsewhere ||Ziaeepour 1997 



Ziaeepour 1999 



Here we use constraint solutions to find a general expression for matching condition on a 
surface of phase transition. This issue has been already discussed in [Peruelle fc Mukhanov 1995 



Deruelle et al. 1997|1 , ||Uzan et al. 199q| . We show that on the light of results obtained 
above, the matching becomes trivial and its physical interpretation more transparent. For 
simplicity, in this section we consider only the case of a flat cosmology. Using formalism 
presented in the previous section, the extension to a curved space-time is straightforward. 
The matter in the early universe has gone through a few number of phase transitions. 
One of the consequences of these transitions is the formation of topological defects which 
could have been the source of initial perturbations (if not completely, at least partially 
|Contaldi et al. 1998| |). To study the evolution of fluctuations after their formation, it 
is usually assumed that transition was very fast i.e. its duration was much shorter than 
evolution time until matter-radiation equilibrium. In this case, a phase transition approx- 
imately defines an equal time 3-space in the space-time. 

The perturbation of the energy-momentum tensor after transition is written as: 



5f, 



(69) 



where ST^ is the perturbation in the ordinary matter (assumed to be relativistic), and 
6^ is the energy-momentum tensor of defects. In stiff approximation, it is assumed that 
Qfj,u evolves separately and is treated as an external source. 

The power spectrum of defects after their formation can only be determined by numerical 
simulations [ [Stephen et al. 1998|| . It needs a complete simulation of quantum processes 
during phase transition and their decoherence which leads to macroscopic classical behav- 
ior of defects. Such a simulation is not yet available and consequently, it is necessary to 
use phenomenological arguments to fix the initial conditions in simulations of perturba- 
tion evolution in presence of defects fVeeraraghavan fc Stebbins 199U[ | , | |Pen et al. 199-1| | . 
Durrer fc Sakellariadou 1997|| . 



It is usually assumed that the effect of defect formation on matter and radiation is tiny 
and consequently perturbative. For studying the evolution of large scale (wavelength) 
perturbations which are important for the formation of large structures today, the lowest 
terms in the expansion of the power spectrum of defects is enough and give an analyti- 
cal expression for initial perturbations [Peruelle et al. 1997j| . Some of amplitudes of the 



expansion terms can be fixed by using the conservation of Mi ,, i.e 0^ = 0. Others are 
usually fixed by some physical arguments, like causality, and/or a matching condition 
that relates a physical quantity in two parts of the space-time separated by the phase 
transition surface. Matching along with other choices define a model for initial perturba- 
tions. These conditions must be at least consistent and treated together. 
For physical reasons like conservation of energy-momentum tensor, the phase transition 
surface is assumed to be a surface of constant density. However, this definition is ambigu- 
ous. In a flat space-time before phase transition, a constant density 3-space is flat. After 
transition, in general, the constant density surface is not equal-time and flat (keeping the 



1 1 



same coordinate definition). To remove this ambiguity, we assume that the isomorphism 
between homogeneous and perturbed space-time manifolds Ml homo and Ai, maps E, the 
transition surface, to a flat equal-time (thus constant density) 3-space. This means that 
at O(0) the density is constant but not at 0(1) (It is the general case. Evidently, it 
is always possible to redefine time parameter such that £ becomes an equal time and 
constant density surface). 

By definition, the space time must be differentiable everywhere including on the initial 3- 
space. This means that a(t), (jij and must be continuous [peruelle fc Mukhanov 1995|[ : 



[^± = and [ify] ± = 0. (70) 

where for any quantity Q, [Q]± = lim e ^ Q(tpr + e ) — Q(tpr ~ e )- Variable t PT is the 
phase transition comoving time. Note that if these conditions are satisfied in one gauge, a 
continuous gauge transformation fulfills them too. Therefore in this sens, they are gauge 
invariant. 

The relation (|70|) is used as the matching condition on the phase transition surface. Here 
we apply the solution of constraints (f|) and (|5p to the matching conditions, and determine 
the unconstrained quantities and K on the initial 3-space. 
From (|10D and ([32]): 



(LW) lj (k, t) = — [Ak 2 ^^ + FT") - 2k c T* c (k i k j + k 2 g ij ) - 

k 2 r)K k 2 n^ 

^(tf* - k -f )]■ Pi) 



and from (0), (0), and (0 



(^horno (^^l/iomo 

(lOSjLo + <LW)l mo + ^W j )^ + g*8K. (72) 



Equation (f4Tf ) gives <50 as a function of 5T*: 



27r^f; + f^(k,t) a 
50= — a' ^ ' 

Dots denote derivative with respect to comoving time t. Before phase transition, Universe 
is homogeneous and the physical metric and extrinsic curvature are equal to their bare 
counterparts. The matching condition for induced metric ( |7D| ) and definition of physical 
and bare metric ( p5|) and (|23|) leads to 50 = and thereafter = g^ = a 2 (t)Sij and: 

5K = -^6f;. (74) 

The matching condition for extrinsic curvature is 5K 1 ^ = 0. In consequence, its three 
independent components i.e. 5K, (LW) 1 ^ (it is 0(1)), and S 1 ^ must separately be zero. 
From (|7D and ([FT]): 

8f* = T* 1 = 0. (75) 



(ff5|) means that there is no total density or velocity perturbation on the initial 3-space. 
If matter is a perfect fluid or a scalar field, i.e. without viscosity, 5Tij = 0. 
The above results show that the matching condition (|70|) is sufficient for determination 
of the initial value of metric and extrinsic curvature and leads to an isocurvature initial 
condition for large wavelengths perturbations for matter without viscosity (e.g. a mixture 
of perfect fluid and a scalar field). No degree of freedom rests to be chosen or fixed by 
other physical arguments. 

Some authors (e.g. ||Pen et al. 1994(1 and [|Turok N. 1996| ) use an ordinarily (in contrast 
to covariantly) conserved pseudo-energy-momentum tensor to define a matching condi- 
tion. In terms of induced metric and extrinsic curvature of the initial 3-space, this tensor 
has the following form: 



5T 00 + £^5K 



5T< 



5f ij - £§ (Hhij - SKy + 8 i:j 8K) 



(76) 



If we impose the continuity of and as before, from ([FJ|) and (|75|) one can triv- 
ially conclude that r^ u conservation condition is fulfilled. In a flat space-time the k ^ 
modes of are zero. [Pen et al. 1994j| and [Turok N. 1996|| use this property as match- 
In contrast to the first prescription, these 
For fixing the rest, Pen et al. 



ing condition and set r 00 = and T 0i = 
conditions don't fix all the independent degrees of freedom 
choose a relation between density perturbation of radiation and dark matter (based on 
the assumption of a white noise power spectrum at superhorizon scales). fTurok N. 1996j| 
chooses these quantities such that at superhorizon scale there is not any perturbation in 
the ratio of different components of the plasma. It is equivalent to an unperturbed initial 
condition. 

In ||Uzan et al. 1998fl , first ([70) is used as matching condition and a series of relations 
between components of metric, extrinsic curvature and energy-momentum tensor of de- 
fects is obtained. Then it is shown that to first order, these conditions is equivalent to 
Too = i.e. the condition used by [[Pen et al. 1994j| and fTurok N. 1996|| for matching. 
From (171), (|75p and ( 75 ) one can immediately and without lengthy demonstration of 
Uzan et al. 1998(1 conclude that continuity of and Kij imposes Too = 0. But the in- 
verse is not true, i.e. if r 00 = and r 0i = 0, it does not necessarily mean that there is no 
total initial perturbation even for a matter without viscosity. Therefore, equivalence of 
two matching prescriptions is one directional and model dependent. Specially, for some 
models like perfect fluid and scalar field, the purely geometric matching of Deruelle et 
al. completely fixes the initial value of the spectrum ( (|74|) and (|75f)), but not the other 
prescription. The advantage of the method presented here is that it is compact and in- 
dependent of the details of the model. This makes the interpretation of the results easier 
and more transparent. 

Finally, we rise the following question: How physically meaningful is matching? A match- 
ing prescription is used in the circumstances that not enough physical information about 
the model is available to fix the initial value of geometry and spectrum. However, even 
a purely geometrical and physically well motivated prescription of Deruelle et al. leads 
automatically to an isocurvature perturbations in a mixture of perfect fluid and scalar 
field. This can be a too much simplification of the reality. One interpretation of this 
result is that any physical process needs a finite time to happen and a complicate process 



1 ■? 



like defect formation can not be replaced by a geometric matching on a 3-space. In this 
situation, it is probably more reasonable to use some phenomenological arguments for 
choosing the initial conditions. 

In conclusion, York method for separation of dependent and independent degrees of free- 
dom in Einstein equations is used to solve analytically the constraint equations for small 
perturbations in space-times with constant curvature. The solution is independent of 
details of the matter model. The results is applied to two commonly used matching 
prescriptions and it is shown that their equivalence is partial and model dependent. 
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